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Abstract 

The  aim  of  this  paper  is  to  outline  the  theory  for  calculating 
the  angular  and  radial  Mathieu  junctions  of  complex 
arguments.  These  functions  are  required  for  the  computation 
of  analytic  solutions  of  electromagnetic  scattering  by  lossy 
dielectric  elliptic  cylinders  and  waveguides.  The 
backscattered  echo  width  of  a lossy  dielectric  elliptic 
cylinder  is  compared  with  the  special  case  of  lossy  circular 
and  weakly  lossy  elliptic  cylinders  and  excellent  agreement 
is  obtained  in  all  cases.  Tabulated  and  plotted  numerical 
results  of  typical  Mathieu  functions  are  presented. 

1.  Introduction 

Many  engineering  and  physics  applications  involve  the 
solution  of  structures  with  elliptical  geometries.  Analytic 
solutions  to  such  structures  require  the  computation  of 
angular  and  radial  Mathieu  functions  in  the  elliptical 
coordinate  system.  Examples  of  such  applications  are  the 
solution  of  elliptical  waveguide  problems  [1-5],  optical 
fibers  [6],  and  elliptical  loaded  horn  and  corrugated 
elliptical  horn  antennas  [2],  field  prediction  inside  lossy 
elliptic  cylinders  for  biological  body  modeling  [7],  and 
scattering  by  mutilayered  dielectric  elliptic  cylinders  [8-12]. 
Furthermore,  analytical  solutions  are  generally  needed  to 
check  the  accuracy  of  numerical  or  approximate  solutions. 
Over  the  years  many  algorithms  and  programs  have  been 
developed  to  address  the  problem  of  computing  Mathieu 
functions  of  integer  and  real  arguments  [4],  [13-22]. 
However,  none  of  these  publications  have  addressed  the 
more  general  problem  of  computing  the  angular  and  radial 
Mathieu  functions  of  complex  arguments.  Caorsi  et.  al. 
have  computed  the  solution  of  electromagnetic  scattering  by 
weakly  lossy  multilayer  elliptic  cylinders  using  first-order 
truncation  of  the  Taylor  expansion  of  each  Mathieu  function 
[23-24].  The  same  authors  have  presented  the  solution  of  the 
scattering  of  an  electromagnetic  wave  by  a lossy  dielectric 
elliptic  cylinder  based  on  the  exact  Fourier  series  solution  in 


terms  of  the  Mathieu  functions  with  complex  argument  [24]. 
Recently,  Amendola  outlined  the  theory  for  the  computation 
of  Mathieu  functions  of  complex  order  by  applying  the 
Floquent  theory  [27].  Unfortunately,  there  are  no 
numerical  results  for  these  functions  given  in  the  literature. 
In  this  paper,  the  solution  presented  in  [15]  is  extended  to 
account  for  the  characteristic  values  (eigenvalues)  of 
Mathieu  functions  of  complex  argument.  Once  the 

characteristic  values  are  computed,  the  complex  Fourier 
coefficients  of  the  Mathieu  functions  can  be  obtained.  These 
coefficients  are  needed  to  compute  the  angular  and  radjal 
Mathieu  functions  of  complex  arguments  and  their 

derivatives  as  will  be  shown. 


2.  Theoretical  Background 


The  homogeneous  wave  equation  (Helmholtz)  in  elliptical 
coordinates  u,  v,  and  z is  given  by 


1 

F2(cosh2  w-cosh2  v) 


3V  dV 

+ 

|iL+*’] 

du2  8v2 

> 

r = o (1) 


The  elliptical  coordinate  system  ( u,v,z ) is  defined  in  terms  of 
the  Cartesian  coordinate  system  (xj>,2)  by 
x = Fcosh(«)cos(v)and  y = F sinh(w) sin(v) , where  F is 
the  semi  focal  length  of  the  elliptical  cross  section. 
Assuming  no  variation  of  the  function  y/  along  the  z-axis. 


the  z-derivative  in  (1)  will  drop  out.  The  resulting 
differential  equation  can  be  separated  into  the  following 
differential  equations  in  terms  of  U(u)  and  V(u).  The 
circumferential  (ordinary)  Mathieu’s  differential  equation  is 
d2V  i dv 2 + (a-2qcos2v)V  = 0 (2) 

and  the  radial  (modified)  Mathieu’s  differential  equation  is 
d2U/du2  ~(a~2q cosh  2u)U  = 0 (3) 

where  <?  = (&f)2  =qR+jql9  k =o)^[ju^€}  and  ex=e\-je\. 

It  can  be  shown  that  there  is  a countable  set  of  characteristic 
values  of  a,  for  which  the  solutions  are  n or  2n  periodic. 
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One  can  easily  distinguish  between  four  kinds  of 
characteristic  values  (eigenvalues)  [15]: 
a = a2r(q) , even  solution,  period  % 

a = a2r+ 1 0?)  9 even  solution,  period  2n 
a = b2r+i  (q)  > °dd  solution,  period  n 
a ~ b2r+\  (4) , odd  solution,  period  2n 
r = 0,1,2,... 

For  every  characteristic  value,  the  solution  of  (2)  can  be 
given  as  a Fourier  series: 

Se2r  (v,  q)=T  A*'  cos  2 mv  (4a) 

m=0 


Selr^  (v,  q)  = £ cos(2 m + l)v 

m= 0 

^t2(v,?)  = |;C  sin(2m  + 2)v 

m-  0 

So2r+t  (v,  q)  = 1 52m+i  sin(2m  + l)v 


(4b) 

(4c) 

(4d) 


Where  the  A’s  and  B’s  are  complex  Fourier  coefficients.  For 
computing  the  Fourier  coefficients,  one  has  to  obtain  first 
the  characteristic  values. 

Substitution  of  the  series  (4a),  (4b),  (4c),  and  (4d)  in  the 
differential  equation  (2),  yields  four  sets  of  equations  for 
computing  the  Fourier  coefficients: 


a2r  ~q  0 
- 2q  anr  - 4 ~q 
0 


0 


0 -q  a2r-(2my 

a2r+\  ~^~q  ~q 


4r 


= o 


(5a) 


-q 

o 


o 

-9  ~q 


0 

- q 

air+\  ~(2m  + \)2 


Afr+] 

A;r+' 


(5b) 


-q  0 
^2r+2“16  ~q 


- (2m  + 2)2 


B]r+2 

n2r+  2 


b'cA 


(5c) 


\.1- 1-9  -?  0 

'5t2r+1‘ 

-9  -<? 

0 

0 

-q 

0 

-q  b2r+x  -(2m  + 2)2_ 

d2»’+1 

. 2m+L 

(5d) 


= o 


We  shall  now  show  the  computation  of  a2r  i in  (5b).  All 
other  characteristic  values  can  be  calculated  in  exactly  the 
same  way.  It  is  clear  from  (5b)  that  a2r-i  can  be  seen  as  the 
eigenvalues  of  the  infinite  tridiagonal  matrix. 


I +q  q 0 
q 9 q 

0 **.  0 
q 

0 q (2m  + 1)2 


(6) 


In  calculating  the  eigenvalues  and  Fourier  coefficients,  it  is 
necessary  to  truncate  the  matrix  after  a finite  number  of 
columns  and  rows  to  obtain  the  required  accuracy. 

Once  the  coefficients  A\rm,  , B^rm++\ , are 
determined  using  the  routine  DEVCCG  from  Microsoft 
IMSL  Math  Libraries,  the  angular  Mathieu  functions  in 
equations  (4)  as  well  as  their  respective  derivatives  can  be 
calculated.  Next,  we  compute  the  radial  Bessel  functions  of 
the  first  and  second  types  using  equation  (3)  along  with  the 
previously  computed  Fourier  coefficients.  The  equations 
used  for  the  computation  of  the  even  and  odd  functions  of 
the  first  type  are 

Jt2r(u,q)  = (-\y4^l2  I(-l Y Alr  (7) 

m- 0 


•W«>?)  = (-l)rVw2 

M _ , 

K-D"  ^2m+l  IVm+1  (*1  Vm  (*2  ) +Jm  (*1  1 (*2  )]  I A,2r+‘ 

m=  0 

Jc2r(u,q)  = (-iyV^72 

Z (- >r Ct4+,(*. )Jm  (x2 (x, )Jm  , (x2 )] / B\r 

m- 1 

Jo2r+](u,q)  = (-])ryfcl2 

Znr  c,  [•/„.,  (*,  vm  ) - j„  o,  ym+,  (*2 )]  / sr 


(8) 

(9) 

(10) 


Jq  Jq 

where  ^ ^—e\x2  = — e’\  By  replacing  Jm  by  Ym, 

we  obtain  the  even  and  odd  radial  functions  of  the  second 
type  as  follows 

Yc2r(u,q)  = {-\y 4^12  s (-1)" 4X (*, V. (*2 VAo  (11) 


Ye2r+](u,q)  = (-\y4^J2 

I (-1)”  C\  [Ym+I  ( * , )jm  (x2 ) +Ym  (x,  )Jm+]  (x2 )]  / A»+l 


(12) 


Yoir(u,q)  = (-iy4^72 

Z (-1)"  *£«,♦,  (X,  )JmA  (x2)  - ,/m+I  (X,  )Km_,  (x2 )]  / B\r 

m= 1 

Yo2r+l(u,q)  = (-iy4^l2 

M 

£ (-  l)m  t YmM  (X,  )Jm  (x2 ) -Jm  (X,  )Jm+i  (x2 )]  / B™ 

m=  0 


(13) 


(14) 


Numerical  results  for  the  derivatives,  with  respect  to  q , of 
the  even  and  odd  angular  and  radial  Mathieu  functions  are 
obtained  using  the  equations  given  in  [14-17]. 
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3.  Electromagnetic  Scattering  by  a Lossy  Dielectric 
Elliptic  Cylinder 

Consider  the  case  of  a linearly  polarized  electromagnetic 
plane  wave  incident  on  a lossy  dielectric  elliptic  cylinder  (of 
permeability  p , permitivity  e and  conductivity  a ) at  an 
angle  fa  measured  with  respect  to  the  X axis  as  shown  in 
[25].  The  electric  field  component  of  the  TM  polarized 
plane  wave  of  amplitude  E0  is  given  by 

(15) 

where  k0  is  the  wave  number  in  free  space.  The  incident 

electric  field  may  be  expressed  in  terms  of  Mathieu 
functions  as  follows 

K = lAexy04)sem{c0,T1)+ 

^ (16) 


Am  =EJm-^-Sem(c0,cos^) 

om  jy  em  (C0  ) om 


Nem{c)=  J [Sem(c,rj)]2dv  (18) 

om  0 om 

where  £ = cosh  u and  t]  = cos  v , c0  =k0F  , F is  the 

semifocal  length  of  the  elliptical  cross  section,  Sem , 
while  Som  are  the  even  and  odd  angular  Mathieu  functions  of 

order  m , respectively,  and  R^  are  the  even  and  odd 
radial  Mathieu  functions  of  the  first  kind,  and  Nem  and 
Nom  are  the  even  and  odd  normalized  functions. 

Electric  field  expansion  inside  and  outside  the  lossy 
dielectric  cylinder  can  be  obtained  by  solving  the  Helmholtz 
equation  in  the  elliptical  coordinate  system  using  the 
separation  of  variables  technique.  The  scattered  electric 
field  outside  the  elliptic  cylinder  for  (^r>^1)  can  be 

expressed  in  terms  of  a sum  series  of  Mathieu  and  modified 
Mathieu  functions  as  follows 

K = X Bem  Rem  (C0  > (C0  > tf)  + 

'”=°  (19) 

iBomyj(c0,^)SoJc  0,»7) 

m=l 

where  Bem  and  Bom  are  the  unknown  scattered  field 
expansion  coefficients,  R™  and  R ^ are  the  even  and  odd 
modified  Mathieu  functions  of  the  fourth  kind.  Similarly, 
the  transmitted  electric  field  inside  the  cylinder  for  £ < 
can  be  written  as 


K = X CeXl(Cvt)Sem(cvl1)+ 

m=0  (20) 

IC0X'yy)Som(C],rj) 

m= 1 

where  cx  =kxF  is  complex,  and  t,  = (0*\nx£x  . 

£x  =€l-j€i  , Cemmd  Com  are  the  unknown  transmitted 
field  expansion  coefficients.  The  magnetic  field  component 
inside  and  outside  the  cylinder  can  be  obtained  using 
Maxwell’s  equations,  i.e. 

«. (21) 

coju}h  8v 

— / dE 

H „ = — (22) 

coju}h  du 

h = F^ cosh2  u - cos2  v (23) 

The  unknown  expansion  coefficients  can  be  obtained  by 
imposing  the  continuity  of  the  tangential  field  components 
at  surface  of  the  cylinder  £ = , i.e. 


I Me,jC]x0)Bem  rf'(c0,6 ) - R^ic^SrMr; 
m=°  Kn(cU l) 


: £ Merm(c^)Aem  R^c0^)^y^-MrR^(c0,4t) 


X Monm(cvc0)Bom  MrR%!\c0,&)- J&>(c U ^7 


= Z MonJciyc0)Aom 

m=0  [ Ro,!(ci’Si) 

Menm  (C0  > C1 ) = lSem  (C\  > V)S„,  (C„  , t])dv 


where  fir  is  the  relative  permeability  of  the  dielectric 
region.  The  TE  scattering  by  a lossy  elliptic  cylinder  is  also 
presented.  In  this  case  a derivation  dual  to  that  of  the  TM 
case  leads  to  the  following  solution  [25] 

I MenJc„c0)B™  erR%\c0,&-  f ? = 


£ Menm(C\,c0)Aem  R^(c^ )S)!Cl’f\)  " £yj(c0,^ ) 


I Monm(ct,c0)B™  erR% (C^)-R^M) = 

r o'  i (28) 

£ Mmn{c^)Aon 

The  scattered  near  and  far  fields  for  the  TM  and  TE  cases 
can  be  calculated  once  the  scattered  fields  expansion 
coefficients  are  known.  The  far  scattered  field  expressions 
can  be  obtained  as  follows 
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El  = 


Le-J*0P  z r [B™  sem  (c0 , n)  + B™Som  (c0 , n)}  (29) 


Kp 


HI  = D-e-**” S r fcX (c0 , 77)  + BllSom (Co , //)]  (30) 
V*0  P 

Far  Field  data  are  usually  expressed  in  terms  of  the 
scattering  cross  section  per  unit  length,  Le.,  the  echo  width. 
For  the  TM  polarization  case  it  is  defined  as 

kf 

r (31) 


am  = 2 Kp  lim 


m 


Eq.  (30)  can  be  put  in  simpler  form  for  the  purpose  of 
computation  as  follows 


v 


[ BmScm  (co , V ) + BomSom  (co , 77)] 

A m 


(32) 


4.  Numerical  Results 

To  determine  the  accuracy  of  the  computer  program,  the 
Fourier  expansion  coefficients  of  Mathieu  functions  are 
computed.  These  results  are  in  excellent  agreement  with 
those  given  in  [16]  for  the  case  of  real  positive  argument  q , 
and  the  maximum  difference  was  found  to  be  1 0*5.  Also,  to 
verify  the  accuracy  of  the  computer  program  for  the  case  of 
complex  argument  q , we  solved  the  problem  of 
electromagnetic  scattering  by  a lossy  dielectric  elliptic 
cylinder  since  no  results  are  available  in  the  literature.  The 
computed  numerical  results  were  compared  with  the  special 
case  of  scattering  by  a lossy  circular  dielectric  cylinder,  and 
with  electromagnetic  scattering  by  weakly  lossy  multilayer 
elliptic  cylinder  [23,26].  Fig.  1 shows  the  numerical  results 
of  the  scattered  normalized  echo  width  a of  a lossy 
dielectric  circular  cylinder  having  a relative  permittivity  of 
€r  =1- y'11.3  and  electrical  radius k0a  = 3.33  for  TM  and 
TE  polarizations.  These  results  compare  very  well  with 
those  presented  in  [26].  In  this  case  the  circular  cylinder 
was  approximated  using  an  axial  ratio  close  to  1,  and  was 
excited  by  an  incident  plane  wave  with  (f>{  = 1 80° . 

Fig.  2 compares  the  back  scattering  normalized  echo  width 
k0<j(0)!4  for  a homogeneous  elliptical  cylinder  with  axial 

ratio  alb-  2 and  £ =4.0.  Two  cases  are  studied;  these 
are  the  lossless  cylinder  case  with  £ = 0.0  and  a weak  loss 

cylinder  case  where  £ = 0.6  . The  calculated  results  are  in 
very  good  agreement  with  those  published  in  [23].  One  can 
easily  recognize  the  reduction  in  the  back  scattering 
normalized  echo  width  due  to  the  presence  of  a small  loss. 
This  can  be  attributed  to  the  fact  that  the  weak  loss  material 
absorbed  part  of  the  incident  wave  energy.  Another  clear 
observation  is  the  shift  in  of  resonance  location.  For  a free 
loss  material,  resonance  is  observed  around  k0a  =3.6,  while 
in  the  weak  loss  case  the  resonance  has  shifted  to  around 


k0a=2A.  Fig.  3 illustrates  the  results  obtained  in  [16]  for 
the  radial  Mathieu  function  of  the  first  kind  and  zero  order 
with  respect  to  u for  different  values  of  q , and  presented 
here  to  validate  the  accuracy  of  the  computer  code.  Finally, 
Figs.  4-13  present  a collection  of  data  for  Mathieu  functions 
(radial  and  angular)  of  complex  arguments.  The  Fourier 
expansion  coefficients  of  Mathieu  functions  for  different 
arguments  are  tabulated  in  Tables  1-4.  These  results  should 
be  very  helpful  for  researchers  solving  problems  in  the 
elliptic  coordinate  system. 

5.  Conclusions 

The  theory  for  calculating  the  angular  and  radial  Mathieu 
functions  of  complex  arguments  has  been  outlined.  The 
computed  backscattered  echo  width  of  a lossy  dielectric 
elliptic  cylinder  was  compared  with  the  special  case  of  lossy 
circular  and  weakly  lossy  elliptic  cylinders  and  the  results 
are  in  complete  agreement.  Selected  numerical  results  for 
Mathieu  functions  are  plotted  and  tabulated  for  limited 
ranges  due  to  limitation  on  space. 
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Figure  1:  TM  and  TE  Normalized  backscattered  echo  width 
versus  <p  of  a lossy  circular  cylinder  with  k0a  = 3.33 and 
£,=1-711.3. 
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Figure  2:  TM  Normalized  backscattered  echo  width  versus 
<p  of  a lossy  dieletric  elliptic  cylinder  with  er  = 4 , a/b^ 

and  ^ =0° 


Figure  3:  Radial  Mathieu  function  of  the  first  kind  and  zero 
order  vs  u for  different  values  of  q. 


Figure  5:  Imaginary  part  of  radial  Mathieu  function  of  the 
first  kind  and  zero  order  vs  u for  different  values  of  q. 


Figure  6:  Real  part  of  the  radial  Mathieu  function  of  the  first 
kind  and  zero  order  vs  u for  different  values  of  q. 


Figure  7:  Imaginary  part  of  the  radial  Mathieu  function  of 

Figure  4:  Real  part  of  the  radial  Mathieu  function  of  the  first  the  first  kind  and  zero  order  vs  u for  different  values  of 

kind  and  zero  order  vs  u for  different  values  of  complex  q.  complex  q. 
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Figure  8:  Real  part  of  the  radial  Mathieu  function  of  the  first 
kind  and  zero  order  vs  u for  different  values  of  q. 


Figure  9:  Imaginary  part  of  the  radial  Mathieu  function  of 
the  first  kind  and  zero  order  vs  u for  different  values  of  q. 


Figure  10:  Real  part  of  even  periodic  Mathieu  function  with 
orders  0-5  and  q = 5 - / 2 vs  v. 


Figure  1 1 : Imaginary  part  of  even  periodic  Mathieu  function 
with  orders  0-5  and  q = 5-j2  vs  v. 


Figure  12:  Real  part  of  odd  periodic  Mathieu  function  with 
orders  1-5  and  q = 5-  j5  vs  v. 


Figure  13:  Imaginary  part  of  odd  periodic  Mathieu  function 
with  orders  1 -5  and  q = 5 - j 5 vs  v. 
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Table  1 : Complex  Fourier  expansion  coefficients  of  even- 
even  Mathieu  function  ( A$m ) 


m\r 

0 

2 

q=5+j0 

Real 

Imag. 

Real 

Imag. 

0 

0.540612446 

0 

0.438737166 

0 

2 

0 

0.65364026 

0 

4 

0.14792709 

0 

-0.426578935 

0 

6 

0 

0.075885673 

0 

q = 5 + j 5 

0 

0.49810151 

0 

0.314442843 

0.171487341 

2 

-0.662750285 

-0.130176082 

0.604942355 

0 

4 

0.168770876 

0.133382581 

0.026159729 

-0.59020582 

6 

flSBBl 

0.089859201 

q = 10  + ylO 

0 

-0.44951258 

0.056260383 

0.029961136 

2 

0.694591697 

0 

0.255064311 

4 

0.746844887 

0 

6 

0.058531831 

0.064924524 

BBSS 

Table  3:  Complex  Fourier  expansion  coefficients  of  odd- 
even  Mathieu  function  ( ) 


m\r 

2 

10 

q = 5+j0 

Real 

Imag. 

Real 

Imag. 

2 

0.933429442 

0 

0.000033444 

0 

4 

-0.354803915 

0 

0.000642976 

0 

6 

0.052963729 

0 

0.010784806 

0 

8 

-0.004295885 

0 

0.13767512 

0 

q = 5 + j5 

2 

0.870294564 

0 

-0.000130329 

m 

-0,396630117 

KBH 

-0.001280155 

0.001222013 

6 

0.028211436 

0.099835145 

BRB 

0.021019319 

8 

0.0053059 

0.132565552 

0.136516812 

<7=10+ yio 

2 

0.734255554 

0 

mn 

-0.000156527 

4 

-0.580164625 

-0.25100546 

HUB 

0.008264692 

6 

0.137051442 

0.19916337 

0.076604876 

8 

0.002015162 

-0.052591626 

0.231256474 

0.260037365 

Table  2:  Complex  Fourier  expansion  coefficients  of  even- 
odd  Mathieu  function  ( ^m+i ) 


j2f+l 

S*2m+1 


m\r 

1 

3 

q=5+J0 

Real 

Imag. 

Real 

Imag. 

1 

0.762463687 

0 

0.077685798 

0 

3 

-0.63159632 

0 

0.30375103 

0 

5 

0.139684806 

0 

0.927728396 

0 

7 

-0.014915596 

0 

-0.201706148 

0 

<?  = 5 + ./5 

1 

-0.517275706 

0.286011205 

-0.041658559 

0.131675807 

3 

0.765785524 

0 

0.239854654 

0.290697818 

5 

■0.197812838 

-0.153451695 

0.878673205 

0 

7 

0.007517743 

0.03807304 

-0.166961766 

-0.194677341 

<7=io+ yio 

1 

-0.433393417 

0.130050198 

0.666818493 

0 

3 

0.775016392 

0 

0.494260893 

■0.299695541 

5 

-0.373034841 

-0.201569455 

-0.01872188 

-0.426395143 

7 

0.054220182 

0.106990887 

-0.184358865 

0.054869005 

Table  4:  Complex  Fourier  expansion  coefficients  of  odd-odd 
Mathieu  function  ( 5^+ \ ) 


m\r 

1 

3 

q = 5+j0 

Real 

Imag. 

Real 

Imag. 

l 

0.940019022 

0 

0.050382462 

0 

3 

-0.336541963 

0 

0.297365513 

0 

5 

0.055477529 

0 

0.931566996 

0 

7 

-0.005089553 

0 

-0.202193639 

0 

q = 5 + j 5 

1 

0.904373155 

0 

0.016113831 

0.092448918 

3 

-0.389707951 

-0.145292103 

0.268912553 

0.288134908 

5 

0.057923929 

0.074976514 

0.87655541 

0 

7 

-0.000091975 

-0.012224892 

-0.16648902 

-0.196783024 

<7  = 10  + j\0 

1 

0.852110863 

0 

-0.043283898 

0.203330605 

3 

-0.460992769 

-0.18202088 

0.234907973 

0.510680316 

5 

0.083408677 

0.142179121 

0.741126923 

0 

7 

0.009774806 

-0.031241645 

-0.081180316 

-0.285889059 
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